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Abstract. Numerical solutions of differential equations are usually not smooth functions. However, they should 
resemble the smoothness of the corresponding real solutions in one way or another. In [2] and [3], a kind of spacial 
smoothness indicators was defined and subsequently applied on the a posteriori error analysis. Here we prove that the 
boundedness of those smoothness indicators is actually a necessary condition for a piecewise polynomial function to 
approximate a smooth function with optimal convergence rate. This should help in validating the error analysis in [2] 
and |3] ■ Moreover, the result of this paper provides an efficient practical method to detect the loss of convergence rate 
due to the lack of numerical smoothness, hence it serves as a criterion for the qualities of many numerical schemes. 
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1 Introduction 

In the literature of numerical solutions of partial differential equations, the smoothness of numerical 
solutions has not been a popular concept. Of course, numerical solutions obtained by using finite 
difference, finite element and finite volume methods are typically not smooth functions. Therefore, 
it seems on the surface that there is no smoothness whatsoever. However, it is widely known that a 
scheme for solving a time dependent problem needs to be dissipative and/or total variation dimin- 
ishing. As a matter of fact, these are actually concerning the smoothness of numerical solutions. 
Although nonsmooth solutions (shocks, interfaces, etc.) are sometimes of great interest, PDE so- 
lutions are usually at least piecewise smooth. A numerical scheme must be able to approximate 
smooth pieces of solutions. Being dissipative and total variation diminishing is certainly necessary. 
However, as shown in this paper, there are actually stronger necessary conditions to be satisfied if 
a scheme is hoped to converge at its desired optimal rate. 

In [2] and [3], a kind of spatial smoothness indicators is proposed. These indicators are sub- 
sequently verified to be bounded in the numerical experiments, which means that the computed 
numerical solutions are "numerically smooth" . Most importantly, these smoothness indicators are 
applied to the local error analysis, playing the role of higher order derivatives. For the scalar 
nonlinear conservation laws studied in [2] and [3], obtaining the local error estimates in terms of 
the numerical smoothness made it possible to do error propagation analysis by directly using the 
i 1 -contraction between entropy solutions. Consequently, a posteriori error estimates of optimal con- 
vergence rates and linear growth were obtained for the RK-DG scheme and the WENO scheme. The 
key advantage of the methodology of [2] and [3] is that the smoothness indicators serve as a bridge 
for bypassing the difficulty of proving any global property of a scheme, caused by the nonlinearities. 



Here, we present the spatial smoothness indicator used in (3J. We prove that the boundedness 
of the smoothness indicator is actually a necessary condition for the numerical solution to converge 
to any smooth function, including the real solution, at the optimal convergence rate. In fact, the 
boundedness of the smoothness indicators is our choice of a global property to deal with in [2] and 
[3j. Since the numerical solutions are computed by the very complex DG, WENO and Runge-Kutta 
schemes, we are certainly unable to give any a priori proof for the boundedness of the smoothness 
indicators. The computation of the smoothness indicators is actually the way to bypass the difficult 
proof. The necessity result of this article confirms that it is reasonable to expect the boundedness 
of the smoothness indicators. The result not only further supports the error analysis of [2] and [3J, 
but also supports the general methodology developed over there. 

Beyond validating the methodology of [5] and [3J, the necessity result provides an extremely 
efficient criterion on numerical solutions of differential equations. Namely, if the smoothness in- 
dicator computed from a numerical solution is too large, then the numerical solution is certainly 
not converging to any smooth function at the desired optimal rate. In another word, whenever 
the smoothness indicators seem too large, either there is some kind of non-smoothness in the PDE 
solution being approximated, or there is something wrong in the numerical scheme being used. 

The result is formulated in a 1-D uniform partition for the simplicity. One can generalize the 
proof to higher dimensions and non-uniform triangulations. However, as the first result of its kind, 
the investigation of a simple 1-D result suffices to show the necessity of numerical smoothness. 



2 The main results 

Let V be the space of polynomials of degree p or less, and Vh be the space of piecewise polynomials 
of degree p or less on the uniform partition a — xq < x\ < ■ ■ ■ < xn = b, h = Xi+\ — Xi = (b — a)/N. 
Let u R be a piecewise polynomial in Vh- The reason for using the superscript R is to indicate that 
u R might be a reconstructed numerical solution, possibly from nodal values, cell averages, or other 
forms of numerical solutions. Define a smoothness indicator for any u R £ Vh as 



with 



and 



S p = (M, D) 

M=(M ,M lr -- ,Mjv-i), where M { = (M°,M/, ■ • • ,Mf), 
D= (£>!,-•• ,D N -i), where A = (£>?, D},--- ,£>?). 



Furthermore, 



and, as in [3J, 



TL/ik d k r, + s k d k R . _. k k k 



D\ = j k /h p+1 - 



In this article, Wj? is the standard notation of the Sobolev space of the functions, where the fc-th 
derivative is L 9 -integrable. H h = W$. See p], Chapter 2. 

Definition 2.1 A piecewise polynomial u R € Vh is numerically W^ 1 -smooth in the partition of 
cell size h, if there is a constant M m such that all the components of S p are bounded by M^; u R 
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is numerically H p+1 -smooth in the partition if there is a constant M 2 such that X^i* ^ [(^i ) 2 + 
(Dj) 2 + • • • + (-Df) 2 ] < M|; is numerically Wf +1 -smooth in the partition if there is a constant 
Mi such that h + \Dj\ + ■ ■ ■ + \D P \] < M x . 

It is obvious that, for the piecewise polynomial u R to approximate a smooth function in the 
interior of all the cells in the domain, one expects \M^\ < 0(1) for all k and i. It is also obvious 
that, if u R is an approximation of a smooth function with the optimal convergence rate, one can 
expect \D®\ < 0(1) for all i. In the case k = p, it is obvious that \D P \ < 0(1) implies that the 
piecewise constant function j^u R (x) has bounded variation. In the main theorem below, we show 
the necessity of the boundedness of the other components of the smoothness indicator. Let's start 
with a simple lemma. 



Lemma 2.2 



Q(D°,D\--- ,D P ) = min 



is a positive definite quadratic form. 



k=Q 



L2(0,i), 



Proof. In the process of calculating the coefficients of the polynomial v = J2k=o Vk£ k to attain 
the minimum, each Vf- must be a linear combination of D° , -D 1 , • • • ,D P , hence the minimum is a 
quadratic form of D°, D 1 , ■ ■ ■ , D p . The quadratic form is positive definite because, if any entry D k 
of (D°, D 1 , ■ ■ ■ , D p ) is non-zero, then we have cither V k + ^ or V k - \^!L ^ o. Due to the 
linear independence of {1, £, • • • , one of the two terms in the minimum must be positive. # 

Theorem 2.3 Suppose that u 6 H p+1 (a, b), and u R is a piecewise polynomial function described as 
above. Then, there is a constant C2 > 0, independent of h, u and u R , such that 



\U-U \\L*(a,b) 



>h p+1 



N-l 



, J2hQ(DlDl,--- ,D p )-C 2 \u\ HP+1(a!b) 
\ »=i 



(2.1) 



Proof. Let V c be the space of piecewise polynomials of degree p or less on the partition a < xi < 
xa < • • • < x N _ 1 < 6, where x { _ 1 = (xi-\ +Xi)/2. Since u G H p+1 (a, 6), there is a constant C2 > 0, 
independent of h and u, and a piecewise polynomial u 1 &V C , such that 



\\U - U 7 || L 2 (a!b) < C2/l P+1 |u|//P+ 1 (a^)- 

Now, since - u R \\ L ^( a ,b) < II" 7 - "IU 2 (a,6) + ||« - u R \\L^( a ,b), we have 

\\u - u R \\ L 2 {a ^ b) > \\u' - u R \\ L 2 {a ^ h) - \\u - ?/||l2 (q ^) 



(2.2) 



> 



\ i=l 



(2.3) 
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,R „„ _ 1 V^P J, 



Let Q* = (Mf + L 4 fc )/2 and w(z) = ELo " Then ' 
x G A+ = and u R - w = —^Y%=a lfr( x ~ x i) k for 1 e ^ — (^i-^^i)- 

= min||w-(u fl -w)||2 2(2; 



mm 



2 
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k=0 



+ 



L2(A7 



1 p J fe 

k=0 



i2 ( A i 



= /i 2p+3 min 



P nfc 



2^ fc! ? 

fc=0 



i 2 ("5,0) 



fc=0 



= ^+ 3 Q(A°,A 1 ,--- un- 
plugging (|2~4| into (f2~3]) . we have proven (|2~T|) . # 



(2.4) 



Theorem 2.4 If u E Wf +1 (a,b), and u R is as described previously, then, there are positi 
stants C\ and C p 2 , independent of h, u and u R , such that 



\\u-u R \\ LHa , b} >h? +1 



\ U \w? +1 (aj 



1<«JV-1 



(2.5) 



Proof. Since u G W p+1,1 (a, b), there is a piecewise polynomial u 7 G V Cl such that 
where C\ > is a constant independent of ft, and it. Obviously, 

II" - U R \\ L i( atV) > Wu 1 - U fl ||jri( 0)6 ) - ||U - U / || i i (a!&) 

:• , v i 



(2.6) 



By the standard scaling argument, and (|2.4[) which remains valid for the current u 1 , it is easy to 
prove that there is a constant Cf 2 such that 



ll«'-« J W< 1; * i+ o >^IK-^IU^ 1 ^)>C^ 3 h^ +2 jQ(DlDl,...,iyt). 

2 2 2 2 v 

Plugging this into (|2.6[) . we have (|2.5p proven. # 

Theorem 2.5 I/u G W^ +1 [a, b], and u R is as described previously, then, there is a constant > 
0, independent of h, u and u , such that 



\\u-u R \\ L ~ {aM >hP +1 



(2.7) 
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Proof. Since u e W r P ) +1 [a, b), there is a piecewise polynomial u 1 e V c , such that 

Hw-^llioo^.b) < Coo/l P+1 |w| W '£+ 1 [a,&]' 

where Co > is a constant independent of h and u. Obviously, 

- 1t R ||i«.( a ,6) > ||u 7 - W fl ||L~(a,f,) - ||« - W i '|U'x>(a,6) 
> max - u fl || ioo( _ 1;X ) - C'o /l ?,+1 |M| W / £ + 1 [a61 



Now 



, let Ui = Wu 1 — u H \\ LoC ( x l X t ), then, by using (|2.4j) for the current u J , 




1 

— (u 1 — u R ) 2 dx 



By plugging into (|ST5]) . we have proven ([27T)1 , # 

3 Conclusion remarks 

According to Theorem 12.31 in order to have 

ll«-« fl ||L»(a,b) <0(F +1 ) 

for any function u £ H' p+1 (a,b), u R must be numerically H p+1 -smooth. That is 

JV=1 

^ hQ(D?,Dl,--- ,Df) < 0(1). 
i=l 

According to Theorem 12.41 in order to have 

\\u-u R \\ LHa>b) <0(h? +1 ) 
for any function u E Wf + (a, 6), must be numerically Wf +1 -smooth. That is 

JV=1 



Y,hjQ{DlD},...,DP)<0{l). 



i=l 



According to Theorem 12.51 in order to have 

\\u-u R \\ L ^ {a , b) <0(hP +1 ) 

for any function u € 14^ +1 [a, 6], u R must be numerically W^ +1 -smooth. That is 

Q(DlDl--- ,Df)<0(l) 



for all i. Because Q(do, d\, ■ ■ ■ , d p ) is a positive definite quadratic form, the last inequality is equiv- 
alent to that \Df | < 0(1) for all k and i. This is what we mean by the necessity of the numerical 
smoothness. To be more explicit, \D%\ < 0(1) is equivalent to \ < 0{hP +1 - k ). That is, the 
jumps of the fc-th derivative of u R need to be as small as 0(h p+1 ~ k ), for k = 0, 1, • • • ,p. 

Although the results and the proofs of the theorems seem to be very simple, the author believes 
that the impact of the theorems could be far reaching. For a numerical solution of a time-dependent 
PDE, the theorem implies that a fully-discrete scheme must enforce numerical smoothness, otherwise 
the convergence of optimal rate must have been lost. Consequently, some of the traditional numerical 
stability notions are possibly inadequate, if they do not enforce numerical smoothness. 
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